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ABSTRACT 

The properties of the spin-down age are investigated. Based on assumption about a uniform 
magnetic field decay law we suggest a new method which allows us to shed light on magnetic 
field decay. This method is applied for following selection: isolated non-millisecond pulsars 
from the ATNF catalog are chosen. Pulsars in the selection are with the spin-down ages from 
4 • 10 4 to 2 • 10 6 years. In order to avoid observational selection we take into account only 
pulsars which are closer to the Sun than 10 kpc. For this selection we restore the uniform 
magnetic field decay law. It appears that the magnetic field decays three times from 4 • 10 4 to 
3.5 • 10 5 years. This function is approximated by modified power-law. We also estimate the 
birthrate of pulsars in our Galaxy and find that it should be about 2.9 pulsars per century. 
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1 INTRODUCTION 

The spin-down ages of radio-pulsars are determined by values of 
P - periods and P - periods derivative. They change with time. If 
n is the braking index and n=3, these characteristics evolve due to 
PP — const law. However, we observe different braking indexes 
in the nature. Therefore we decide to investigate a uniform function 



f(t) — v PP /olBq. This function is not a constant because of 
magnetic field decay, alignment or some combination of both these 
phenomena. Although we suppose that fit) changes because of 
magnetic field decay, this assumption does not affect generality of 
our results. 

Indeed, there are a few evidences both theoretical and obser- 
vational which indicate that magnetic field of pulsars should decay. 
The theoretical evidence is that the current keeping the magnetic 
field should decrease because of following reasons: finite conduc- 
tivity of matter in neutron s tars (ohmic dissipation), the Ha ll ef- 
fect or ambipolar diffusion (iGoldreich & Reiseneggerlll992l) . The 
observational evidence is that temperature of n eutron stars is con- 
nected with its magnetic field jPons et al.ll2007h . M oreover, kinetic 
age d eterminations also show magnetic field decay dHarrison et al.l 
1 19931) . Last studies find that ages determined by expansion of the 
parent supernova nebula are often smalle r than those determi ned by 
the spin-down of the associated pulsar dZhang & Xidl201ll) . Also 
pulsars are ofter close r to the galactic plane than it can be expected 
( Guseinov et al ] |2004l). Though the population synthesis still leave 



problem unsolved jPopov et al.ll20ld : [Faucher-Giguere e t al. 2006; 
llgoshev & KholtvgirJl201ll) . 

Meantime, there is an additional effect of the magnetic field 
decay which still has not been carefully studied. In fact, changes 
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in magnetic field affect the spin-down age of radio pulsar. Many 
authors define the spin-down age as: 



P_ 

2P 



(1) 



This formula estimates the age of the pulsar based on assumption 
that one steadily increases its period P from to current value. Al- 
though this quantity is commonly used, we understand its proper- 
ties poorly. Indeed, this estimation should be the same as the real 
age only for quite old pulsars (with ages more than 10 4 years). 
Moreover it strictly depends on magnetic field decay and align- 
ment. For instance, if B, magnetic field strength, decreases follow- 
ing a power low with exponent 7, the spin-down age increases as 
r = yt (see the second Section for deductions). 

Nevertheless, B is exchanged to r does not give us additional 
information. We still need to observe r(t) for every particular pul- 
sar to restore magnetic field history. This difficulty can be avoided 
if we assume a uniform magnetic field decay law B = Bof(t) 
where Bo is the initial magnetic field and f(t) - is universal time- 
dependent law. This assumption leads us to conclusion that all pul- 
sars can be sorted by r. It intends that the higher r is, the higher 
real age is also. Therefore we can consider pulsar with particular 
age r as a stage in the life of every pulsar in the Galaxy. 

After all, it is appeared that r is B-dependent and vice verse. 
So, does it intend that we lost only available age estimation? Prob- 
ably not. Let us have a look to different process which like annual 
rings in cross-section cut of a trunk deliver us age. It is quantity 
of pulsars. This quantity is determined by a birthrate, conditions 
which leads to pulsars detection and life-span of them. The first of 
them is number of pulsars which on average are born during fixed 
period of time. It is quite constant value which hardly ever vari- 
ate significantly due to fixed life-span of massive stars and absence 
of starbursts in our Galaxy. Two other values strongly depends on 
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pulsars characteristics. However, on average for ensemble it is seen 
equilibrium. Then let us continue analogy with annual rings. Ev- 
ery one thousand years new pulsars are born as a new annual ring 
during summer. When we count quantity of pulsars with some par- 
ticular spin-down age we notice lacks and excesses which mark 
changes in magnetic field. 

Consequently, there are two independent age estimate r and 
taver (deduced from birthrate of pulsars). The spin-down age r is 
affected by changes in B. Nevertheless, many readers can suggest 
a few reasons why r is bigger than t even without magnetic field 
decay. For instance, we observe smaller quantity of pulsars with 
bigger r because they are weaker. Although this effect is signifi- 
cant, it can be avoided if we take into consideration a part of our 
Galaxy. Also effect of alignment can increase the spin-down age. 
Both of these two effects will be discussed in detail. 

The article is structured as follows. The second Section is de- 
voted to clarification of some important features of r. In the third 
Section we discuss uniform law of magnetic field decay and con- 
sequences of this assumption. The fourth Section is devoted to a 
derivation of kinetic equation which allows us to restore real ages 
of pulsars. In the fifth and seventh Sections we present discussion 
of key selection effects. The method and results are presented in the 
sixth Section. 



(a cos 2 f + b sin 2 £)d B 2 = (a + (b - a) sin 2 O^B 2 = PP (3) 

This expression means constant braking rate with weak angular de- 
pendence (\b — a | < a). For sake of simplicity elsewhere in this 
article except the seventh section we use a = b. 

Meanwhile no angular dependence involves no alignment. 
The alignment is a phenomenon when £ decreases with time. 
Therefore taking into consideration expression (0 in form J3]) we 
assume that there is no alignment or its effect is negligible on inter- 
val of time shorter than 10 6 years. 

The expression (0 is not so easy-to-use as integrated expres- 
sion for purposes of our article. We can treat (0 as a differential 
equation, and its solution is: 



P 2 (t) 



I a{r')B 2 {T')dT' + C 
Jo 



(4) 



Where C can be found from boundary condition P(0) = Po where 
Po is the initial period of the pulsar. Let us combine ([T} and ©: 



T(t) = 



f*a(T , )B 2 (T , )dr , + 0.5P 2 



(5) 



a(t)B 2 (t) 

This is an integral expression for the spin-down age which helps us 
to study its properties. 



2 MATHEMATICAL PROPERTIES OF r 
2.1 Connection with magnetic field 

Let us write the standard expression for magnetic field if a pulsar is 
braked due to electromagnetic radiation: 



aB z = PP 



(2) 



Where a is a constant for the given pulsar at given moment 
which depends on mass, momentum of inertia and angle £ be- 
tw een magnetic pole and the axis of rotation. In early work 
bv lOstriker & Gunnl dl969l) the authors considered pulsar brak- 
ing as a consequence of dipole radiation emission in vacuum. 
However, iGurevich et al.l dl993l) showed that the pulsar has to 
posses an extend ed magnetosphere fille d with plasma. There- 
fore mecha nism dOstriker & Gunnl [l969|) does not work in it- 
self. In factlGureyich et al.l d 19931) deduced similar formula © as 
lOstriker & Gunnl (1 19691) did, except the angular dependence. The 
former believes that a ~ cos 2 (£) when the latter supposes that 
a ~ sin 2 (£). 

On the other hand astronomers hardly ever measure both 
B and £ in their observation. The angle £ i s measured by be- 
havior of polarization angle during a pulse dManchesterl Il97ll : 
iLvne et al.lll97ll : iMcCulloch et al.lll978l : iMorris et al.lll98ll) . The 
direct magnetic field measurements are possible by detection of 
synchrotron absorption features in high energy part of spectra 
(iBlandford & Scharlemann Il976l) . However, this measurements 
tend to differ from quantity (|2j at some constant factor. It is might 
connected with different heights where high energy and radio emis- 
sion photons are bor n. The height is smaller, m agnetic field is 
stronger. In the paper dTauris & Manchesterlll998l) it is shown that 
there is a relation although it does not behave like sin 2 £ or cos 2 £, 
it is much more weaker. Therefore the conclusion about angular 
dependence can not be deduced from observations right now. 

Meantime, in observations we see that all isolated pulsars are 
braking. Therefore we suppose angular dependence in the follow- 
ing form: 



2.2 The simplest magnetic decay models 

Currently astronomers tend to rely on a power-law and exponential 
magnetic field decay models. The first describes behavior of the 
field in a form: 

t V 



»-*(£)" 

Using we find the spin-down age behavior: 
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1 



p2 



(6) 



(V) 



T(t) ° + 2aB(t) 2 7 + l t+ 2aB 2 (t) 

Clearly seen that the spin-down age in this model is proportional 
to real age and the rate is constant. More often astronomers use 
the exponen t ial ma gnetic field decay law in pulsar population e.g. 
Pop ov etallbOlGh . Let us suppose that the field decays as follow- 
ing: 



B = Boe~ t/T 
Then we apply relation (|5]l: 

r (t) \- — 

K) e- 2 */r ^2aB 2 (t) 



P§ = T __ T _ 2t/r | Po 
2m 2 2 e 



2aB 2 (t) 



(8) 



(9) 



Here the spin-down age is growing exponentially. These two the 
simplest examples hint that the spin-down age is highly sensitive 
to behavior of B and is weakly affected by the initial period and 
initial magnetic field. It explains why we choose this quantity in 
our study. 

2.3 Averaging 

About 1700 isolated pulsars are known in our Galaxy. Each of 
them has specific parameters set. This set contains magnetic field 
strength, period, period derivative, radio-luminosity, rotation axis 
orientation etc. Some of the characteristics are connected with each 
other. For instance, formula (0 connects period, derivative of pe- 
riod and magnetic field. Using this formula we can deduce any of 
this characteristics by others. 
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There are characteristics which are not connected with each 
other. They are independent. For instance, nowadays we can not 
suggest how initial magnetic field of a neutron star and initial pe- 
riod are related. Therefore let us call two values independent further 
if we can not suggest any physical law which associates them and 
we do not see any correlation in observational data. 

The ensemble of pulsars contains about 1700 discrete objects. 
Every object there has the whole set of observed characteristic. 
Thereby, we observe the discrete set of characteristics. Meantime, 
there is no physical argument in favor that it is not possible for 
an additional pulsar to have some intermediate characteristics. It 
means that observed ensemble is a realization of some distribution. 
Therefore it is quite natural to replace one by continuous distribu- 
tion and work further with it instead of particular pulsars. 

Let us suppose that £ is a property of a pulsar such as a pe- 
riod, magnetic field or something like this. Then we introduce the 
distribution function oj((). Let us choose it as a normalized one: 



J a 



j(C)d( = 



(10) 



When we average some quantity on this distribution we get an ex- 
pectation value of this quantity. Due to presence of variety in obser- 
vation data for every parameters, we use expectation value in our 
analysis. 

2.3.1 Averaging on initial magnetic fields distribution 

As it is mentioned above the initial magnetic fields distribution and 
initial periods distribution are independent. Therefore we can aver- 
age on these characteristics independently. Let us consider Bo as £• 
We also introduce here the following designation: 



rB 2 

r{t)\ Bo = / r(t,B' )uj(B f )dB 
Then the expression ([5} can be rewritten: 



r(t)\ Bo 



J^a(r f )B 2 (r f )dr f + 0.5P§ 



a(t)B 2 {t) 



(11) 



(12) 



The initial magnetic field might affects a only during first about 
10 3 years when it significantly deforms n e utron stars if the field 
is strong enough ([T hompson et al.l |200(J; iG hosh & Chakrabartv 
201 L 2009; Haskell 2008; IColaiuda et al.ll2008l : IWassermanll2003l : 
Ostriker & Gunnlll969l) . On the contrary the aim of our investiga- 
tion are middle-aged pulsars. Therefore we can carry a out of the 
integral. Moreover as any physical quantity B(t) and oj(Bo) are 
continuous, so it is possible to exchange integrals: 

B 2 



r(t)\ Bo = 



B 2 (t\ Bp) 
B 2 (t,B ) 



j(B )dBodT f + 



2aB 2 



(13) 



Here we take assumption about behavior of magnetic field. We 
suppose that B(t) — Bof(t), and f(t) is a monotony function. 
Thereby: 

2 f 4j7^^) dB ^' + -= (I 4 ) 

Pit) 2<xB*P(t) 



r(t)\i 



nt rB 2 
JO J Bl 



Because f(t) does not depend on the initial magnetic field and us- 
ing relation dTTTb we take: 



SU 2 iT')dT> 
Pit) 



2aBlP{t) 



(15) 



It is again the spin-down age with a small disturbance. This dis- 
turbance will be estimated after the averaging on the initial periods 
has been done. 



2.3.2 Averaging on initial periods distribution 

Let us do here the averaging on initial periods distribution simi- 
larly as it is done on initial magnetic fields in the previous section. 
Firstly, consider ( as Pq. We designate the averaging by: 
rP2 

r(t)\p = / r(t,P Q f )oj(P f )dP f (16) 
J Pi 

Using the relation $5$ we can rewrite: 



T(t)\ Po 



J t a(T')B 2 {T')dT' + 0.5Pi 
a(t)B 2 (t) 



(17) 



The initial period might affects on a only during first about 10 3 
years when it significantly deforms neutro n crust if it is very short 
dOstriker & Gunnlll969l : ICutler et al.ll2003h . after this time gravita- 
tional radiation brakes neutron stars enough that changes in a can 
be neglected. As it is mentioned above the aim of our investigation 
are middle-aged pulsars. Therefore we can carry a out of integral 
due to a is constant on our time-scale. Moreover, as any physical 
quantity, P(t) and cj(Po) are continuous, so it is possible to ex- 
change integrals: 

^ = [ C ^W^ dP ^' + 2^ (18) 

As it is mentioned above initial magnetic field and initial period 
are not connected. Therefore we can carry out J Q aB 2 {r') from 



T(t)\p - ^oTT\ H 



(19) 



P{t) 2aB*f>(t) 
And again we get similar the spin-down age with small disturbance. 

2.3.3 Averaging on both distributions 

Now, when we already know expressions ([T5t and (|T9), we average 
on both characteristics simultaneously: 



r(t)\ 



Pol Bo 



r(t)\ Bo \ 



The result is: 



r(t)\p , Bo 



P(i) 



p2 
^0 



2aBlP{t) 



(20) 



(21) 



We should carry f(t) out of parenthesizes because it affects sim- 
ilar on the first and the second summand. Then we take typ- 
ical values for members in numerator and denominator from 
dFaucher-Giguere et al.lE006h as P = 0.1 sec, 5 = 4- 10 12 G 
and a = 3.2 • 10 -19 . Therefore this disturbance member can be 
estimated as: 



10" 



2.5 • 10 11 sec ; 



10 3 years (22) 



p2 

aBl = 10- 39 • 4 • 10 25 

It means that t\b ,p ~ r without any disturbance members for 
quite old pulsars. Consequently our method can be applied for pul- 
sars with the spin-down age which is larger that 4 • 10 4 years. This 
value defines boundary where our method can be applied. 
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3 THE UNIFORM MAGNETIC FIELD DECAY LAW 

In order to clarify how physical the assumption about the uni- 
form law is we need to consider what properties of a neutron star 
can affect the magnetic field decay. These properties are rotation, 
crust conductivity, accretion, temperature and properties of mag- 
netic field itself such as configuration and strength. Let us discuss 
these factors. 

If we exclude objects which accrete matter, we consider only 
pulsars with similar properties. Indeed, majority of neutron stars are 
born with short periods with small variance (IFaucher-Giguere et al.l 
2006). Neutron stars are born in supernovae explosions, and con- 
sequently should have similar characteristics in beginning because 
of similar conditions of supernova explosion. We exclude accretion 
since we take into consideration only isolated pulsars with quite 
long periods. As for magnetic field configuration, pulsars tend to 
brake such way that dipole magnetic field configuration is looked as 
the most probable. Nethertheless, pulsars with the strongest mag- 
netic fields brake in a different way. They emit radiation due to 
destruction of magnetic field. These pulsars are magnetars. How- 
ever, we exclude them from our selection either. Therefore this as- 
sumption about the uniform magnetic field decay law seems quite 
natural. 

Moreover, this assumption also gives us time order. Indeed, 
if we have two pulsars one with the initial magnetic field B\ and 
another one with P2, then for every time t\ > tiT\ > T2 and vice 
verse. Indeed: 

^aBjf(t')dt' £fvw , i;;f(t')dt' (23) 



Tl = 



aBlf(U 



f(ti) < f(t2) (f(t) is monotonous function) it means that the first 
summand is bigger or the same as T2 . Simultaneously the second 
summand is positive, and consequently n > T2 

When we study the ensemble of real pulsars we should re- 
member that there is a condition when a pulsar turns off. This con- 
dition is determined by combination of B and P. It means that pul- 
sars with lower fields turn off earlier. Due to this fact, our selection 
is not complete for larger ages. 

Currently the death-line for isolated radio-pulsar s is deter- 



Currently the death-line tor isolated radio-pulsars 
mined by following condition ( Bhattacha rva et al .11 1992b : 



B 

p2 



< c 



(24) 



Here c is a constant. On the other hand we can rewrite this expres 
sion via B and r: 
1 



r(t)> 



2aB cf(t) 



(25) 



Then we estimate B = 10 13 G, c = 1.7 • 10 n G s~ 2 , f(t) = 0.1 
and get r tU rn-off = 10 8 years. However, near 10 6 — 10 7 years ef- 
fects connected with alignment begin influencing, see I Young et al.1 
(2010) and references therein. 



4 KINETIC EQUATION AND STATISTICAL AGE 
ESTIMATION 

For the whole pulsar ensemble we can write a kinetic equa- 
tion. Let us consider small bits at the spin-down age and 
real age space. In this two dimensional space we con- 
sider real ages as time-coordinate and the spin-down ages 
as space-c oordinate. T hi s app r oach differs from th a t sug - 
gested by iBeskin et all dl986h : IPhinnev & Blandfordl dl98lh : 



IVivekanand & Naravanl Jl98lh : lDeshpande et al.l (fl995) where au- 
thors considered distribution function out of (P, P, £). In our work 
we suppose that there is no apparent angular dependence, the mag- 
netic field B changes due to uniform law B = Bof(t). Therefore r 
is a function only of t for quite old ages (see the second Section for 
deductions). In these terms n(r, t) is pulsars distribution function. 
It denote quantity of pulsars with the spin-down from r to r + dr 
and real age from t to t + dt. We can write now the continuity 
equation for the whole pulsar ensemble: 



dn ( dr 

— — V div ri- 
ch V dt 



dr^ 



U-V 



(26) 



Here U is a source of pulsars and V is a summand which describes 
pulsars death. The simplest source of pulsars in this coordinate 
space should be ribirthrate^(O) which sets the spin-down age for 
all new-born pulsars as zero. In reality, majority of pulsars are born 
with the spin-down age near 10 3 years. And the total birth rate is 
determined as: 



R 



p(B )p(Po)dB dP 



(27) 



Here p(B) determines the initial distribution of pulsars by their 
magnetic fields, and p(Po) by their initial periods. We can define 
parametrization Bo = Po(A,r) and Po = Po(A,r) which de- 
termines this way that VA P 2 /(2aPo) = r. Then we integrate 
p(Bo) p(Po) by this parameter A: 



P(r) 



R 



P 

Jo 



p(B (\))p{P (\))d\ 



p(t)cIt 



(28) 



(29) 



It is important to mention that there is the highest possible 
value of r at the birth. And pulsars with r > 4 • 10 4 years at the 
birth are rare. Therefore we can truncate this distribution function 
at value 4 • 10 4 years. Because we take into consideration pulsars 
with the spin-down ages smaller than 2-10 6 years we can neglect V. 
We also suppose that the pulsars ensemble is stationary. Following 
these assumption we can rewrite (|26t : 



dn dr , N dB 



(30) 



In order to find n(r, t) in a certain position we should integrate this 
relationship: 



dt 



The result is: 



l(r, t) = nbirthr. 



dt_ 



(31) 



(32) 



It is important to remind that n is density, so the quantity of pulsars 
with the spin-down age smaller than r is: 



(33) 



N(r,t) = / n(r,t)dr 
Jo 

Let us put (ITTt into d32l and integrate the result: 

N(t) 



t(r) = 



^birthrate 



(34) 
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5 SELECTION EFFECTS 

Results of the method, either it is right or not, might be spoiled 
due to significant observational selection. In fact, we observe only a 
small part of all galactic pulsars. And we do not have any solid base 
to suppose that this small part correctly reproduces characteristics 
of the whole galactic ensemble. Indeed, selection effects are able to 
hide some group of pulsars from the whole ensemble. In this case 
selection effects change an average and a variance of the observed 
distribution. 

For purposes of our study the most important selection is that 
which is caused by decrease of radio luminosity for old pulsars. 
Astrophysics suppose that radio luminosity of pulsars is a constant 
part of their rotational energy. Due to this hypothesis old pulsars 
with longer period should have smaller supply of kinetic energy, 
and consequently they should be weaker radio sources. We are able 
to observe we aker sources only clo ser to the Sun. Moreover, fol- 
lowing article dGurevich et al.l l 1993) pulsars with higher magnetic 
fields brakes quicker, and consequently have smaller luminosity. It 
might cause apparent magnetic field decay even if it does not have 
place. Our method should consider full selections for all spin-down 
ages before pulsars shutdown. Therefore this effect demands care- 
ful consideration. 

On the other hand it is difficult to estimate this selection on the 
theoretical ground because of our poor understanding of the radia- 
tion mechanism. Therefore we draw our attention to observational 
data. We choose isolated radio-pulsars which are not mill i secon d 
ones from the ATNF pulsar catalogue Q Ma nchester et al.l ([2005). 
We divide them into four groups: (the spin-down age is in years) 
r e [10 3 — 2 • 10 5 ], t £ [2 • 10 5 — 5 • 10 5 ], r £ [5 • 10 5 - 2 • 10 6 ] and 
r £ [4 • 10 7 — 10 8 ] and plot them in FigureQ]due to their distance 
from the Sun. It is clearly seen that shares first three of them are 
similar when we consider distances which smaller than 10 kpc and 
the last shows evidences of observational selection. Consequently, 
this procedure lets us to avoid observation selection. 



6 DESCRIPTION OF THE METHOD 

The first step of our method is to plot n(r) as a cumulative distri- 
bution. This distribution is presented in the Figure [2] When we plot 
this figure we extract from the ATNF pulsar database only isolated 
non-millisecond pulsars which lie closer than 10 kpc to the Sun. 
It allows us to avoid observational selection (see the fifth Section 
for details) and effects of accretion. Therefore both our assump- 
tions about the uniform magnetic field decay law and fullness of 
selection seem correct. Then we exchange discrete distribution to 
continuous one by approximation of polynomial, and treat this dis- 
tribution as t\p ,b = t the universal spin-down distribution (see 
the second Section for proof). We choose polynomial of 6th degree 
to avoid too detailed, wave-like behavior near 10 3 — 10 4 years and 
on the other hand to carefully reproduce the form of distribution 
from 10 4 - 2 • 10 6 years. We fit distribution by last- square method 
and get following formula: 



g(x) — wx 6 + ax 5 + bx 4 + cx 3 + dx 2 + fx + g 



(35) 



With coefficients a = 0.00883, b = -0.14894, c = 1.2942, d 
-6.2086, / = 16.54, g = -21.665 and w = -2.102 • 10" 4 



http ://www. atnf . csiro . au/research/pulsar/psrcat/ 




Figure 2. Cumulative pulsar distribution by the spin-down ages. Solid line 
- distribution of pulsars and dashed line is approximation. 



The second step of our method is conversion from n(r) to 
t(r). This can be done accordingly to the fourth part of our arti- 
cle. This conversion needs some estimate of pulsar birth-rate. We 
do this estimate rely on consequences of our model field decay. We 
suppose that f(t) is a monotony function. Therefore for every mo- 
ment of time r > t. We set r = t at the moment 4 • 10 4 years and 
calculate birthrate due to this assumption. Then we relay on fact 
that it is possible to observe only 10% of all pulsars and comparing 
volume of 10 kpc Sun-centered circle to the whole volume of the 
Galaxy we recalculate birthrate from 1 . 1 visible pulsars per thou- 
sand of years to 2.9 pulsars per century. It looks like a minimum 
estimate because in reality it appears that r(t) ^> t. However, some 
astronomers believe that the magnetic field decay occurs on time- 
scale much more than 10 4 years. It means that when age of pulsar is 
about 10 4 years his real and the spin-down ages are similar. There- 
fore this estimate should give the real birthrate of pulsars in the 
Galaxy. Following this method we get nb irthrate — 2.9 pulsars 
per ce ntury which is slightly greater than iFaucher-Giguere et al.1 
(2006) got. This birth-rate is two times greater than 1.4 pulsars 
per century from IVranesevic et al.l (120041) estimate. Nethertheless 
our estimation als o agrees with the rate of supernova explosions 
(iKeane & Kramedl 20081) . 

Let us write here the equation (|2H : 



T(t): 



f 2 (t) 



(36) 



Let us carry f(t) from denominator to left part and differentiate by 
t: 



V{t)r(t) + f(t) • (f(t) - 1) = 
Here it is possible to separate variables: 

/(*) _ m , i 

fit) 2r(t) + 2r(t) 
It can be integrated: 



fit) 



(37) 



(38) 



(39) 



Accordingly to this relationship we numerically calculated values 
of f(t) for every moment of time. The result is presented in Figure 
[3] It is possible to approximate this curve by following expression: 
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Figure 1. Distributions of pulsars by distance for different intervals of the spin-down ages a) 10 3 — 2 • 10 5 years, b) 2 • 10 5 — 5 • 10 5 years, c) 5 • 10 5 — 2 • 10 6 
years, d) 4 • 10 7 - 10 8 years 



Here 7 = 1.17, a = 0.034, c = 0.84, to = 10000 years. 
7 DISCUSSION OF ALIGNMENT 

Alignment is a phenomenon when the angle between the magnetic 
field pole and the rotational axis decreases with time. It is extremely 
important for our work because strong alignment affects on the 
spin-down age. It even can ruins our conclusion about absolute- 
ness of r. It happens because pulsars with high r might be young 
pulsars with £ — > 90°. Therefore we may mistakenly mix young 
pulsars with older ones. Then we can not say that the spin-down 
ages distribution some way reflects magnetic field decay 

In order to protect our method from the alignment support- 
ers criticism we need to consider alignment as a part of modern 
braking models. One of them is model of magneto-dipole radiation 
dOstriker & Gunnlll969l) . In this model: 



240000 440000 640000 840000 1.04e+06 1.24e+06 




150000 200000 
Real age, years 



PP = aB 2 (t) sm 2 x 



(41) 



Figure 3. The magnetic field decay law f(t) (solid line) for middle-aged 
radio pulsar with approximation (dashed line). 



There is an invariant which remains constant during angle evolution 
jEliseeva et al.ll2006h : 



Spin-down age: the key to magnetic field decay 1 




5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 
Real age, years 

Figure 4. The spin-down pulsars evolution plotted by their real age. Differ- 
ent lines refer to different initial angles. Solid line - the initial angle is 85°, 
dashed - 80°, big dots - 75°, small dots - 70° 



cos x 



(42) 



Using commonly known relationship sin 2 £ + cos 2 £ = 1 we can 



rewrite (I4TL (l42l into differential equation: 

P(£)P(t) + I 2 dp P(tfaB 2 (t) = aB 2 (t) 



(43) 



The numerical solution of this equation lets us understand the 
value of alignment. We calculated the spin-down ages due to equa- 
tion fl3t with typical initial magnetic field (B = 5 • 10 12 G) 
and initial period (Po = 0.1 sec). As we can see from the Fig- 
ure |4]only pulsars with the initial angles greater than 75° at 40000 
years have t he spin-down age near t his value. On the other hand, 
this work bv lOstriker & Gunnl (11969 ) supposes homogenous angle 
distribution. Therefore only (90 — 75) /90 = 16.7 per cent of all 
pulsars appears in the area where we measure the birth-rate. How- 
e yer, our birthrate estimat e is in good agreement with other works. 
If Ostriker & Gunnl |l969) had been right, we would have taken six 
times smaller birthrate estimate. Therefore we believe that there is 
no alignment or its effect is negligible. 



8 CONCLUSION 

Based on the assumption about the uniformity of magnetic field 
decay law we suggest a new method which allows us to restore 
the behavior of magnetic fields on the surface of middle- age ra- 
dio pulsars. We prove that this method can be applied to analysis 
of statistical properties of r. We restore the magnetic field decay 
law and find that it can be fitted by modified power-law d40t . The 
birthrate of pulsars in the Galaxy is estimated based on the assump- 
tion that r « t at 40000 years. This estimate is in good agreement 
with other ones. Moreover, r(t) evolution is calculated taking into 
account alignment. It is found that if alignment had taken place, 
we would have taken at once three times smaller birthrate by our 
method than it is. 
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